www.gusucode.com > 基于matlab编程蒙特卡罗仿真设计源码程序 > 基于matlab编程蒙特卡罗仿真设计源码程序/第五讲 蒙特卡洛方法的应用/litiR2.m
function Rguji=litiR2(t,thetaa1,thetaa2,thetab1,thetab2,mm) %t 是要求统存寿 %thetaa1 是元件A1的数学期要 %thetaa2 是元件A2的数学期要 %thetab1 是元件B1的数学期要 %thetab2 是元件B2的数学期要 %mm 是随机实验次数 frq=0; randnuma1 = exprnd(thetaa1,1,mm); randnuma2 = exprnd(thetaa2,1,mm); randnumb1 = exprnd(thetab1,1,mm); randnumb2 = exprnd(thetab2,1,mm); for ii=1:mm if (randnuma1(1,ii)>t)|(randnumb1(1,ii)>t) pass1=1; else pass1=0; end if (randnuma2(1,ii)>t)|(randnumb2(1,ii)>t) pass2=1; else pass2=0; end if (pass1*pass2)==1 frq=frq+1; end end Rguji=frq/mm